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Abstract 

We report the first calculation of the simplest but most fundamental transverse momentum 
dependent (TMD) distribution of quarks in the nucleon, i.e. the time-reversal-even unpolarized 
TMD quark and antiquark distribution with isoscalar combination, within the framework of the 



i-S^ ' chiral quark soliton model. The nonperturbative account of the deformed Dirac-sea quarks within 

pi • the theoretical scheme enables us to make a reliable predictions not only for the quark distribution 

, but also for the antiquark distribution. We found that the predicted average transverse momentum 

^ square (/c^) of quarks and antiquarks depends strongly on their longitudinal momentum fraction 

\^ • X, which means that the frequently used assumption of factorization in x and k± is significantly 
OO ! 

QQ . violated. It is also found, somewhat unexpectedly, that the average transverse momentum square 



of antiquarks is considerably larger than that of quarks. 



I. INTRODUCTION 



Growing attention has recently been paid to the transverse-momentum-dependent (TMD) 
parton distributions also called the unintegrated parton distributions [l^. The reason of 
particular interest stems from the fact that they play important roles in the theoretical 
description of single spin asymmetries in various hardprocesses including semi-inclusive 
deep inelastic scatterings, Drell-Yan processes, etc. [sl-j^l. The TMD parton distributions 
are interesting themselves also from more general viewpoint. They are expected to give a 
three-dimensional view of the parton distributions in momentum space, thereby providing us 
with complementary information besides what can be obtained through generalized parton 
distributions [sl -j^. 

It has been established that at the leading twist there are totally eight TMD parton 
distributions - Among those, six are even under naive time-reversal transformations 
(for brevity, we call it T-even), while the remaining two are T-odd 12|- ij]. To meet the 
general requirement of factorization and gauge-invariance, these TMD parton distributions 
must be defined with the gauge link operators called Wilson lines 15| - 20| . Since these 
Wilson lines, which can be interpreted as simulating the effects of gluon initial and final 
state interactions, are process dependent IGj, the TMD parton distributions are in general 



non- universal 



2l|, 



22l |. Still, a remarkable difference between the T-even and T-odd TMD 
parton distributions should not be overlooked. For the former, the presence of the Wilson 
lines is less crucial in the sense that they do survive even without such gauge links. On the 
contrary, for the T-odd TMD parton distributions, the introduction of the Wilson lines is 
of fatal importance since they do vanish identically in the absence of such effects of initial 
and final state interactions The existence of the former types of TMD distributions 

appears only natural since the intrinsic motion of quarks described by the naive bound state 
wave function without any complex phase is likely to be three dimensional and it should be 
described by some TMD distribution function. 

Unlike the integrated parton distributions, our knowledge on the TMD parton distribu- 
tions is very poor. Very recently, through the analysis of azimuthal asymmetries both in 
hadron-hadron collisions and in semi-inclusive deep-inelastic scatterings, Anselmino et al. 
extracted a T-odd TMD distribution called the Sivers function [23]. There also exists 
an attempt to extract another T-odd distribution called the Boer- Mulders function I2J]. In 
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these studies, however, a flavor independent Gaussian distribution of the transverse momen- 
tum under the factorized ansatz in the variables x and k± is assumed by reason of simphcity, 
although there is no compelling reason for this choice. Since enough empirical information 
to test such assumption cannot be expected at the present moment, any information from 
models of baryons and/or the lattice QCD would be extremely valuable. 

Although there exist a lot of model calculations for the integrated parton distributions, 
there are not so many for the TMD distributions. (See 25|] for a compact overview of 
those studies.) The leading- twist T-even functions were calculated in a spectator model 
with scalar and axial- vector quarks 26|], and in a light-cone quark model where the Fock 
expansion is truncated to consider only 3 valence quarks . On the other hand, the T-odd 



functions are investigated in the spectator model with scalar diquarks 
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291] and with 



34| . and in the 



scalar and axial- vector diquarks 301 ]. 3l| . 25| in the MIT bag model 32| - 
constituent quark model |35|. Very recently, the flrst lattice QCD calculation of the lowest 
moments of the TMD quark distributions fi{x,kj_) and giT{x,kj_) has also been reported 
361 ] . (Although our main interest here is concerned with the TMD quark distributions inside 
the nucleon, we recall that there also exist some investigations on the TMD distributions in 
;he pseudoscalar mesons including the pion based on the SU(3) Nambu-Jona-Lasinio model 



37| -|39i.) 

Now, the purpose of the present paper is to evaluate the simplest but most fundamental 
TMD parton distributions in the nucleon, i.e. the T-even unpolarized TMD quark and 
antiquark distributions with isoscalar combination, within the framework of the chiral quark 
soliton model (CQSM) [4^,j4l|. The greatest advantage of the CQSM over the other models 
of baryons is its full account of nonperturbative chiral dynamics of QCD, which has been 



proven to be essential for the physics of light quark distribution functions [42| - 48| . The 
nonperturbative account of the deformed Dirac-sea quarks within the scheme enables us to 
make a reliable predictions not only for the quark distributions but also for the antiquark 
distributions, as already been proven through its good reproduction of the famous NMC 
(New Muon Collaboration) observation, i.e. the dominance of the d sea over the u one in 



the proton, etc. 



491 ]. 461 ]. 501 ]. We therefore expect that the same model would provide us 



with valuable and reliable information also on the nature of the transverse motion of quarks 
and antiquarks inside the nucleon. 

The plan of the paper is as follows. First, in sect. II, we derive theoretical formulas 
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necessary for the calculation of the TMD quark and antiquark distributions within the 
framework of the CQSM. Next, in sect. Ill, the predictions of the CQSM are shown and the 
detailed discussion on them will be given. Conclusion of our analysis is then given in sect. IV. 



II. UNPOLARIZED TMD QUARK DISTRIBUTIONS 

We are interested here in the simplest but most fundamental TMD parton distribution 
function, i.e. the unpolarized TMD quark distribution. The unpolarized TMD quark distri- 
bution function of flavor a in the nucleon, averaged over its spin, is defined as 

where 

^0,^] = Ve-'3l,^ M^yd^^ (2) 

is the so-called gauge link operator, also called Wilson line, connecting the two different 
space-time points and ^, by all possible ordered path. This gauge link is known to simulate 
the interaction of the outgoing quark field with the spectators inside the hadron, and it works 
to ensure the color gauge invariance of the above definition of the TMD quark distribution 

Here, we want to evaluate the TMD quark distribution within the framework of the 
CQSM. Since the CQSM is an effective quark theory, which does not contain gluonic degrees 
of freedom at least explicitly, we can drop the gauge link operator. As a matter of course, 
this simple procedure is fatal for predicting the T-odd TMD distributions like the Sivers 



function 12|, 13| or the Boer-Mulders function [14], since the complex phase arising from the 
gauge link is indispensable for the existence of such T-odd distribution functions. Roughly 
speaking, the T-even distribution function we are to calculate below can be thought of as 
a "static" TMD quark distribution, which can be computed directly from the bound state 
wave functions of a target nucleon without accompanying complex phase arising from the 



effects of final state interactions 



5l|. 



Though the definition ([T]) of the TMD distribution function itself is Lorentz-frame inde- 
pendent, it is convenient to evaluate it in the nucleon rest frame, where P = and Pq = Mn 
with Mjv being the nucleon mass. This enables us to rewrite the exponential factor e^^'^ 
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contained in ([T]) in the following manner : 



(3) 



, Just similarly to the integrated unpo- 
42I - 44 1, 46], 48 1 , the isoscalar distribution 



Here, we have used the relation = — (or ^"^ = 0), together with the standard definition 
of the light-cone vector, = (,^° ± S,'^) / 

To proceed further, we notice that the large- iV^ behavior of the unpolarized TMD quark 
distribution depends on the isospin combination 53] . 
larized distributions extensively studied before 
is dominant in the large-A^^c limit, and survives at the mean-field level, while the isovector 
combination appears as the first order correction in the collective angular velocity operator 
f2 of the rotating soliton. (We recall that is an XjNc quantity.) Since the evaluation of 
the 0(^2^) contribution is much harder, let us concentrate on the calculation of the isoscalar 
part in the present paper. Using the formalism developed in the previous studies 42] - 48] , 
we find that the isoscalar unpolarized TMD quark distribution /"^'^(x, fc^) in the CQSM is 
given in the following form : 



J J (2vr)^ 



n&occ 



(4) 



with = 7° 7^. Here, $„(a;) is the eigenf unctions of the Dirac hamiltonian H with the 
corresponding eigenvalue En, i.e. 



with 



H 



OL ■ V 



(5) 



(6) 



In (jlj), sunin^occ denotes the summation over the occupied states in the hedgehog mean 
field. In deriving (j3]), the projection into a nucleon state with given momenta P has been 
achieved, as usual, by integrating over all shift of the center-of-mass coordinate R of the 
soliton, 

(P'] •■■ ]P) = J rf3i2e*(^'-^)-^ ■■■ . (7) 
Introducing the eigenfunctions in the momentum representation, 



, I px 



(2vr) 



*n(p), 



(8) 



can readily be transformed into the form : 

to J2 



E ^'^"''^ I ~^n^P) (1 + ~^n{p). (9) 



X 

n(^occ 



Now, by noting that e*^'^ = e*^ ^ e*^^'^^ = e *p ? e*^-^'^^, one can carry out the 
integration over and to obtain 

ndocc 

As usual, the numerical evaluation of the above expression is carried out by using the dis- 



cretized momentum basis of Kahana and Ripka 



the smeared distribution defined by [43| 



53|,[5J], so that it is convenient to introduce 



f:;^\x,k^) = \ dx'e ^ r+\x,k^). (n) 

After some straightforward algebra, we therefore get 

/OO 1 f E +k \ ^ 

dhe'^y^'- ) <i.|^(fcx,A:3)(l + a=^)<l„(fc^,A;3). (12) 

Here, we have changed the integration variable from to k^. We recall now that is 
simultaneous eigenf unctions of the hamiltonian if, and the grand spin operator K = J + 
and its projection Mx on the z-axis. The functions can be expanded in terms of the 
discretized plane- wave basis of Kahana and Ripka as 



= E E ^^"^ (13) 

where the index a distinguishes 4 independent plane-wave basis with definite grand spin 
Mk, and parity. (We point out that the expansion coefficients c"^^ can be chosen real 
numbers.) On the other hand, i labels the discretized momenta ki, which is determined by 
imposing the boundary condition 

jKihD) = 0, (14) 

at the radius r = D chosen to be sufficiently larger than the typical soliton size. (The 
number of momentum bases is made finite by introducing the maximum momentum kmax 



such that ki < kmax-) Since the Kahana-Ripka plane- wave basis is given in the spherical 
representation, we also express k± and in the spherical coordinates as 

(k^)^ = ki sin 9 coscp, (15) 
{k±)y = ki sin 9 sin0, (16) 
ks = ki cos 9. (17) 

Since the final answer for the TMD distribution f^~^'^{x, k±) should be independent of the 
azimuthal angle of the transverse momentum k^, the following replacement is justified 

cPl{l + as)cPpj — I ^0L(l + «3)</>ft-. (18) 
The result of integration can schematically be written as 

^ (1 + "3) </>/3i = 5ij J^apiki, COS 9), (19) 

where J-'aj3{ki, cos 9) is a function of ki and cos6'. The Kronecker delta 6ij arises from the 
fact that the operator (1 + 03) does not change momentum. Because the magnitude of the 
transverse momentum k± = \k^\ is given externally, the relations f|T6l) - f|T7j) dictates that, 
once the plane- wave-basis momentum ki is given, both of k^ and cos 9 are fixed as 

fca = ±^kf-kl = ks,^, cos9 = ^ = cos9,. (20) 

We therefore approximate the expression fill I) as 

f^^\x,k^) = N^J2J2I1 ^„,(fc.,cose,)-^e-^(^-^)'. (21) 

n(iocc a, (3 i v T 

In the next section, we show that this algorithm in fact works if we take large enough bound- 
ary radius D, at least for the calculation of Dirac sea contributions. Unfortunately, we find 
that the same algorithm does not work very well for the calculation of the discrete valence 
level contribution (corresponding to n = 0). This causes no serious problem, however, since 
we can use simpler method to evaluate this part. We just go back to the expression ffTOj) 
before introducing smeared distribution, and obtain 

f:y{x,k^) = MNN,^ik^,k3 = xMN-Eo)il + a3)%ik^,k3 = xMN-Eo), (22) 

where $0 and Eq are the eigenfunction and the eigenenergy corresponding the valence level. 
Since $0 is a discrete bound state wave function anyhow, we can evaluate the above expres- 



sion without any difficulty by using its momentum space wave function as used in 



45| 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The basic lagrangian of the model contains two physical parameters, the weak pion decay 
constant the dynamically generated effective quark mass M. As usual, /,r is fixed to be 
its physical value, i.e. f-,^ = 93 MeV. On the other hand, M is taken to be 375 MeV, which 



is favored from our previous analysis of the nucleon spin structure functions [47j,[48|. 

The model contains ultraviolet divergences so that it must be regularized by introducing 
some physical cutoff. Following the previous studies, we simply use the Pauli-Villars regu- 
larization scheme with single subtraction. In this scheme, any nucleon observables including 



quark distribution functions in the nucleon are regularized through the subtraction 



(O)-^ = {Of - [-^\of^\ (23) 

Here {Of denotes the nucleon matrix element of an operator O evaluated with the original 
effective action with the mass parameter M, while (O)^^^ stands for the corresponding 
matrix element obtained from {Of by replacing the parameter M with the Pauli-Villars 
cutoff mass Mpy. Demanding that the regularized action reproduces the correct normaliza- 
tion of pion kinetic term in the corresponding bosonized action, Mpy is uniquely fixed by 
the relation 

J^M'log^ = U. (24) 

For M = 375 MeV, this gives Mpy — 562 MeV, leaving no adjustable parameter. 

A short comment may be necessary for the regularization scheme explained above. As was 



shown in 55|, the Pauli-Villars scheme with a single subtraction term is not a completely 
satisfactory regularization procedure. It fails to remove ultraviolet divergences of some 
special quantities like the vacuum quark condensate, which contains quadratic divergence 
instead of logarithmic one. For obtaining finite answers also for these special observables, the 
single-subtraction Pauli-Villars scheme is not enough. It was shown that more sophisticated 



Pauli-Villars scheme with two subtraction terms meets this requirement [55|]. Fortunately, 
the self-consistent solution of the CQSM obtained in this double-subtraction Pauli-Villars 
scheme is only slightly different from that of the naive single-subtraction scheme, except 



when dealing with some special quantities containing quadratic divergences [55|]. Consider- 
ing the fact that the calculation of quark distribution functions, much more the TMD quark 
distribution functions, in the CQSM, is extremely time-consuming and that the most nucleon 



observables are rather insensitive to which regularization scheme is chosen, we shall simply 
use here the single-subtraction Pauli-Villars scheme. (The use of more time-consuming dou- 
ble subtraction scheme is mandatory, however, for some special parton distribution functions 
containing quadratic divergences. An example is the chiral-odd twist-3 unpolarized distri- 
bution functions e(x) investigated in [56| -|59|.) 



P+''{x, yfcx) with rc = - 0.1 P^'^ix, k^;) with a; = + 0.1 




— bare data 

— smooth fit 



0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

i-2 1.2 

FIG. 1: The Dirac-sea contributions to the unpolarized TMD distribution with isoscalar combina- 
tion f'^^'^{x, kj_) for two typical values of x, i.e. x = — 0.1 (left panel), and x = + 0.1 (right panel). 
In both figures, the dashed curves represent the bare numerical predictions of the model obtained 
with D = AO/ M, k 

max — 10 M , and 7 — 0.3, while the solid curves are their smooth fits. Note 
that the dashed and solid curves are almost indistiguishable in the right panel. 

First, let us check whether the numerical algorithm proposed in the previous section 
in fact works. As pointed out there, there is no difficulty in the calculation of discrete 
valence level contribution to the TMD quark distribution, so that we concentrate on the 
Dirac-sea contributions. Shown in Fig{T] are the predictions of the CQSM for the Dirac-sea 
contributions to /"^"'(x, k±) for two typical values of longitudinal momentum fraction, i.e. 
X = — 0.1 (left panel) and x = +0.1 (right panel). In these figures, the dashed curves are 
the bare numerical predictions of the model obtained with D = 40 / M, k^ax = 10 M, and 
7 = 0.03. The observed fluctuating behavior, which is a little stronger for the negative value 
of X, is due to the use of the discretized basis with fairly small 7. We fit these results with 
smooth function shown by solid curves. 
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After demonstrating that our numerical method in fact works as long as one takes large 
enough boundary radius D, let us now show the full CQSM predictions for the unpolarized 
TMD quark and antiquark distribution functions. We recall that the function fc^) 
with positive x can literally be interpreted as quark distribution, whereas the function with 
negative x should be interpreted as antiquark distributions with an extra minus sign as 

r+^(-x, k^) = - r+\x, k^), (0 < X < 1). (25) 

Figj2]show the CQSM predictions for f'^^'^{x, k±) as functions of for 6 typical values of 
X, i.e. X = 1.0 X 10-^ 0.1, 0.2, 0.4, 0.6 and 0.8. On the other hand, Figj3] shows fc^) 
for 6 different values of x, i.e. x = — 0.8, — 0.6, — 0.4, — 0.2, — 0.1 and — 1.0 x 10"^. In these 
figures, the contributions of the discrete valence-level quarks and those of the Dirac-sea 
quarks are illustrated by the dashed and dash-dotted curves, respectively, while their sums 
are denoted by the solid curves. First, let us look into the quark distribution with x > 0. 
As expected, the contributions of Dirac-sea quark is most important in the lower x region, 
while the contribution from the discrete valence level dominates over that from the Dirac-sea 
quarks as x increases to approach 1. Somewhat unexpectedly, the contributions from the 
Dirac-sea quarks are seen to have longer range tail in k± space, as most clearly seen from 
the figure corresponding to x = 0.1. (The plot of /""'"'^(x, A;_l) is not the best way to see this 
unique feature. A better way to see it more clearly would be to examine the graph of the 
/i;_L-weighted distribution /""'""'(x, kj_), the integral of which over kj_ gives the integrated 
quark distribution /"+'^(x). See the discussion later, for more detail.) 
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FIG. 2: The CQSM predictions for the unpolarized TMD distribution with isoscalar combination 
in the positive x region. They are shown for 6 typical values of x, i.e. x = 1.0 x 
10^^,0.1,0.2,0.4,0.6 and 0.8. In these figures, the dashed and dash-dotted curves respectively 
stand for the contributions of the discrete valence level and of the deformed Dirac-sea quarks, 
while their sum is represented by the solid curves. 
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FIG. 3: The CQSM predictions for the unpolarized TMD distribution with isoscalar combination 
in the negative x region, which correspond to the antiquark distributions, according 
to the rule /""'"'^(x, k±) = — f^~^'^{— x, k±) with < x < 1. They are shown for 6 typical values of 
X, i.e. —0.8, —0.6, —0.4, —0.2, —0.1 and —1.0 x 10^^. The meaning of the curves is the same as in 
FigEl 
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Turning next to the distribution with x < 0, corresponding to the antiquark distribution, 
we find that the contributions from the Dirac-sea quarks dominate over that of valence 
level. Here, the fact that the Dirac-sea contributions have higher k±_ components than the 
valence-level quark is much more evidently seen. An importance notice here is that the 
total contributions shown by the solid curves are all negative for any x. In consideration of 
the extra minus sign indicated in (!25|) . this means that the positivity of the TMD antiquark 
distribution /""'"'^(x, k±) is legitimately fulfilled in the CQSM. We emphasize that the proper 
inclusion of the Dirac-sea contribution is crucial for the fulfillment of this fundamental 
property of unpolarized parton distribution. It is this unique feature of the CQSM that 
enables us to make a reasonable predictions not only for the quark distribution but also the 
antiquark (or the sea quark) distributions. 



r+''{x,kA_) : CQSM /fc_L) : Gaussian 




FIG. 4: The contour plot for the isoscalar unpolarized TMD quark distribution function 
/"+'^(x, /cj_). The left panel corresponds to the prediction of the CQSM, while the right panel 
to the schematic distribution in the factorized form /"^'^(x, fc_|_) = /"+'^(x) ^ ^ e~ ^'^ ^ ^^'^'^ with 
(kl) = 0.25 GeV^. 

For lack of enough empirical information, the TMD quark distribution is often assumed 
in a factorized form with the Gaussian distribution in k± as 

r+'(x,fcx) = r^\x) ^-j-e-'l/^'l\ (26) 

71 {k^} 
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with {k]_) ~ 0.25 GeV^ Then, it may be of some interest to compare the predictions of the 
CQSM for the TMD quark distribution /""'"'^(x, with this factorized ansatz. We show 
in FigJH the contour plot for the isoscalar unpolarized TMD quark distribution function 
/"^'^(x, fcx). The left panel corresponds to the prediction of the CQSM, while the right 
panel to the factorized form given by fl2B]) . (For the integrated distribution f^~^'^{x) in the 
factorized form, we use the prediction of the CQSM, i.e. it is obtained from the CQSM 
prediction for f'^~^'^{x, k±) after integration over k±.) One clearly sees significant difference 
between the behaviors of the two distributions especially in the lower x region, x < 0.4, 
where the magnitude of the TMD distribution is dominantly large. 



/"+'*(a;, /fc_L) : CQSM p+^{x,k_L) : Gaussian 




FIG. 5: The contour plot for the isoscalar unpolarized TMD antiquark distribution f^~^'^{x,k±). 
The meaning of the curves is the same as in FigUl 

Figj5] show a similar comparison for antiquark distribution f'^~^'^{x, k±). One again ob- 
serves significant difference between the prediction of the CQSM and the simple parametriza- 
tion (!26l) . According to the prediction of the CQSM, the TMD antiquark distribution with 
lower X is seen to extend over higher k± region, as compared with the schematic form. At 
any rate, the analysis above indicates that the frequently used assumption of factorization in 
X and kj_ is most probably violated. Later, this statement can be made a little more quanti- 
tative through the investigation of the average transverse momentum square of quarks and 
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antiquarks in dependence of the longitudinal momentum fraction x. But, before doing it, 
here demonstrate realistic nature of the predictions of the CQSM, through the analysis of 
the empirically well-known integrated quark and antiquark distributions, the information of 
which is hidden in the TMD distribution /"^'^(x, k±). 
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FIG. 6: The integrated unpolarized quark and antiquark distribution function f'^~^'^{x) with 
isoscalar combination obtained from the corresponding TMD distribution via (j27p (left panel). 
The meaning of the curves is the same as in Fig.l. Shown in the right panel is the unpolarized 
quark and antiquark distribution with isovector combination, which has been evaluated directly, 
i.e. without reference to the TMD distribution, in the previous papers. 



To demonstrate the reliability of the CQSM predictions for the unpolarized TMD quark 
distributions shown in Fig.l and Fig. 2, we evaluate from f"^~^'^{x, k±) the corresponding 
integrated quark distributions through the relation. 



d'k^r+\x,k^). 



(27) 



The answer is shown in the left panel of Figj6l which nicely reproduces the previous results 
for u{x) +d{x) calculated directly without referring to the TMD quark distribution [48]. (A 
finer look however reveals that our new result for f'^^'^{x) shows a slight enhancement in the 
magnitude of the Dirac-sea contribution as compared with our old result shown in Fig. 1(a) 



of 



481], for which the check of the model space dependence, i.e. the D ^ oo and k^, 



oo 



limit, was not satisfactory enough. The new result for the d{x)/u{x) ratio to be shown 
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below is a little sensitive to this change of the isoscalar unpolarized parton distribution.) 
One confirms again that the proper introduction of the Dirac-sea contribution is essential for 
reproducing the positivity condition of the unpolarized antiquark distribution. One also sees 
that, in the lower x region, even the quark distribution is dominated by the contributions of 
Dirac-sea quarks. If one compares this figure of f^^'^{x) with the Fig|2]for TMD distribution 
fc_L) especially at x = 1.0 x 10~^ and x = 0.1, one can confirm that the longer 
range tail in k± of the Dirac-sea contribution is an important factor leading to its strong 
enhancement in the small x region, which we observe for the integrated distribution /""'"'^(x). 
In the right panel of Figl6], we also show the previously calculated isovector unpolarized quark 
distribution /""'^(x), since it is necessary for the following comparison of the d{x)/u{x) ratio. 
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FIG. 7: The CQSM predictions for the difference function d{x) — u{x) in the proton in comparison 



with the FNAL 



6l| and HERMES data [t^] (left panel), and for the ratio d{x)/u(x) in comparison 



6l| and NA51 data 



60] (right panel). 



with the FNAL (Fermi National Accelerator Laboratory) 
The curves labeled CTEQ4M, MRS(R2) and MRST are phenomenological PDF fits to the data 
before and after the E866 measurement. 



In order to demonstrate realistic nature of the CQSM predictions for the antiquark distri- 
butions, we show in FigHits predictions for the difference u{x)—d{x) and the ratio d( x)/ u(x) 



in comparison with the corresponding empirical data. As in the previous works 



47j,|48j, the 
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theoretical predictions for the distribution functions at the high energy scales are obtained 
by solving the DGLAP equations at the next-to-leading oder, with the predictions of the 
CQSM as initial scale distributions given at the scale (5^„j = 0.30 GeV^. (Here, the suffix 
"ini" should be understood as an abbreviation of the word "initial".) A good agreement 
observed in the left panel means that the model reproduces the famous NMC observation, 
i.e. the dominance of the d-sea over the u-sea inside the proton, without any artificial fine- 
tuning. Shown in in the right panel of FigI7]is the CQSM prediction for the ratio d{x)/u{x) 
at = 30 GeV^ in comparison with the old NA51 data [60] and the newer E866 data ex- 
tracted from the neutrino scatterings. The new MRST parton distribution function (PDF) 



fit including the E866 data j6ll together with the older CTEQ4M and MRS(R2) PDF fits 
are also shown for reference. One sees that the prediction of the CQSM is qualitatively 
consistent with the E866 data as well as the new MRST fit. We emphasize again that it is 
a completely parameter free prediction of the model. 

Now, after demonstrating that the CQSM well describes the fundamental physics of 
unpolarized sea quark distributions, we come back to its predictions on the unpolarized 
TMD distributions for quarks and antiquarks. From the already-given unpolarized TMD 
distributions /"^'^(a;, k±) as functions of x and k"^, it is straightforward to evaluate the 
average transverse momentum of quarks and antiquarks as a function of x : 

_ jd^k,kir-%x,k) 
^"^^""'^ ~ jd^ku-+%x,k) ■ ^^^^ 

The resultant {^^^{x)) is shown by filled circles in FiglHl The solid curve here is a smooth 
fit to the numerical results by an Sth-order polynomial as 

8 

(fci(x)) = 5^c„a;'^ (29) 



n=0 



where 



Co = 0.311122, ci = -0.536064, cs = -2.4806, 
C3 = 18.7331, C4 = -41.7892, cs = 40.5805, 

C6 = -11.3483, C7 = -7.48546, cg = 4.26044, (30) 
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for X > 0, while 



Co = 0.311122, ci = -0.805606, Ca = 14.9032, 

C3 = 114.34, C4 = 376.97, C5 = 693.994, 

ce = 733.782, C7 = 415.798, Cg = 97.7558, (31) 

for X < 0. 




FIG. 8: The CQSM prediction for the average transverse momentum square of quarks (x > 0) and 
antiquarks (x < 0) as a function of the longitudinal momentum fraction x. 

One clearly sees that the average transverse momentum square is strongly dependent on 
the longitudinal momentum fraction x of quarks and antiquarks. This reconforms that the 
frequently used assumption of factorization in the variables x and k± is significantly broken. 
(The plausible breakdown of the usual Gaussian ansatz for a factorized k± dependence was 



also indicated in the model calculation by Pasquinni et al. 27|]) Very recently, the first 



lattice QCD st udy of the transverse momentum distributions of quarks in the nucleon has 



been reported 361]. The validity of factorization in x and k± was examined, and it has been 
concluded that the factorization hypothesis holds within the statistics of the simulation in 
sharp contrast to our present analysis based on the CQSM. In our opinion, the results of the 
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lattice QCD simulation at the present level cannot be taken too seriously by several reasons. 
First, the simulation is preliminary in that it was performed in the fictitious heavy-pion 
region around ^ 500 MeV. What is lacking here is the nonperturbative chiral dynamics 
of light quarks, which is the core of the celebrated NMC observation, or the physics of 
Gottfried sum. Second, the test of factorization is far from persuasive, since the evaluation 
of full x-dependence is beyond the scope of the presently-known algorithms of lattice QCD. 

Turning back to the x-dependence of the average transverse momentum square of quarks 
and antiquarks, we continue to analyze the prediction of the CQSM. Let us first look into 
the positive x region corresponding to the quark distribution. Very strangely, the average 
transverse momentum square has a minimum around x ~ 0.25, and it increases as x becomes 
larger. This feature can be understood as follows. Around x ~ (0.2 — 0.3), the integrated 
quark distribution f^~^'^{x) has a peak, which is dominated by the contribution of the valence- 
level quarks. This means that the relative importance of Dirac-sea quarks, which has higher 
k± components, as compared to the valence-level quark, is smallest around there, which 
explains the above-mentioned feature of (fc^(x)) as a function of x. One may also notice 
that the average transverse momentum square around x = is fairly large, which indicates 
that the quarks (and antiquarks) with small longitudinal momentum fraction would carry 
sizable amount of orbital angular momentum along the z-axis, i.e. along the the direction 



of nucleon spin. We recall that this observation is consistent with the results of [6^, in 
which a direct calculation of quark orbital-angular-momentum distributions was carried 
out. (See also the analyses of the quark orbital angular momentum from the viewpoint of 
the generalized parton distribution functions [63| -[66l|.) 

Another unique feature of the CQSM prediction is that the magnitude {k']_{x)) is much 
larger in the negative x region, corresponding to the antiquarks. Another way of convincing 
this feature is to compare the two quantities defined below : 

{klf ^ lo dxjklix)) r^%x) ^ ^^2^ 



/o 

•<0 



Kdxf-+^ 



X] 



_ j\ dx{kl{x))r+\x) _ jl dx {kl{x)) f-^\x) 



-'^ix) 



j\ dx f-+^{x) /; dx p 

which represent the average transverse momentum square for quarks and antiquarks, respec 
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tively. Numerically, we find that 

{klf = 0.224 GeV^ (34) 
{klf = 0.445 GeV^, (35) 

which clearly shows that the average transverse momentum of antiquarks is much larger 
than that of quarks. We can also estimate the average transverse momentum square of 
quarks and antiquarks altogether from 

(A;i)«+« = -^-^ \ )'J Li (36) 

^ ' J\ dx f^+'^ix) ^ ' 

which gives 

{kl)Q+Q = 0.266 GeV^ (37) 

This value of average transverse momentum square is remarkably close to that used in 



the phenomenological analysis of the semi-inclusive reactions 23|]. Note, however, that the 
average transverse momentum is a scale dependent quantity, which is believed to grow as 



increases. The predictions of the CQSM corresponds to the scale ~ (0.30 — 0.40) GeV^ 



47l | , [48| , while the phenomenological analysis in [23|] should correspond to somewhat higher 
energy scale. In any phenomenological analysis, one must keep in mind this scale dependent 
nature of the average transverse momentum square of quarks and antiquarks. 

Before ending this section, we make a short comment on the so-called transverse- 
coordinate representation of the TMD parton distributions. The Q^-evolution of the TMD 
parton distributions or the unintegrated parton distributions is described by the Catani- 



Ciafaloni-Fiorani-Marchesini (CCFM) equations [67|]-[69|. A nice feature of the CCFM equa- 
tion at the leading order is that they reproduce the convectional leading-order Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations for the integrated distributions. The 
CCFM equations are known to take a particularly simple structure if one introduces the 
so-called transverse-coordinate representation of the parton distributions [t^ -73|. To 
explain it, we represent below the TMD quark (or antiquark) distribution of flavor a as 
f^{x,k±,Q'^) with k± = \k±\, where the dependence on the scale is also shown explic- 
itly. The transverse-coordinate representation /'*(x, b, Q"^) of a quark (or an antiquark) of 
flavor q is introduced as the two-dimensional Fourier transform of the corresponding TMD 
distribution f"-{x,k^,Q^) : 

r{x,b,Q') ^ I d'k^e^'-"- r{x,k^,Q'), (38) 
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where b is the transverse coordinate vector. Clearly, at 6 = 0, the function f°'{x,h,Q'^) 
reduces to the familiar integrated distribution f"-{x,Q^). The greatest advantage of using 
the transverse-coordinate representation is that the corresponding evolution equations for 
the unintegrated distributions, f]\fs{x,b,Q^),fs{x,b,Q'^) and fg{x,b,Q^), become diagonal 
in the transverse coordinate b. Here, just for illustration, we show the evolution equation 



for the nonsinglet (NS) distribution. It takes the form 7^, 75| : 



Q'^ fNs{x,b,Q') = / dzP,,{z) (39) 







x[Q{z-x)M{l-z)Qb]fMs[^,b,Q') - fNs{x,b,Q')], [40) 

which certainly is diagonal in the variable b. Once these evolution equations are solved in 
a similar way as the standard DGLAP equations, the average transverse momentum square 
for a given distribution at any can be evaluated from 



,2^^a _ I d'k^klnx,k^,Q') _ ^ ^nx,b,Q' 



(41) 

6=0 



r{x,b^) = ——e-''^-^^H-{x,^ = 0,t = -A^). (42) 



The transverse-coordinate representation of the unintegrated parton distribution discussed 
above should not be confused with more popular impact-parameter dependent parton dis- 
tributions introduced by Burkardt , although there should be strong correlation between 
them. The simplest impact-parameter dependent parton distribution function /"(x, b^) is 
the two dimensional Fourier transform of the unpolarized generalized parton distribution 
function H"'{x,C,,t) with zero skewdness parameter ^ = : 

We recall that the frequently used parametrization [76| of H°-{x,^ = 0,t) is of non 
factorizable form in x and t as 

H^ix,^ = 0,t) = r(x)exp[/3'^(x)t], (43) 

with 

p%x) = [a' Hl/x) + Ba] {1-xf + Aax{l-xf, (44) 

which means that /°(x, b±) is also non-factorizable in x and b±. This in turn strongly indi- 
cates that the transverse-coordinate representation of the unintegrated parton distributions, 
and consequently, the TMD parton distributions are most likely to be non-factorizable in 
the variables x and k±, as our predictions based on the CQSM shows explicitly. 
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IV. CONCLUSION 



To conclude, we have reported the first calculation of the simplest but most fundamental 
TMD parton distribution in the nucleon, i.e. the unpolarized TMD quark and antiquark 
distributions with isoscalar combination, within the framework of the CQSM. The reahstic 
nature of the CQSM predictions was demonstrated by showing that the integrated unpolar- 
ized distribution /"^'^(a;), obtained from the corresponding TMD distribution f^~^'^{x,k±), 
is consistent with the available empirical information on the light-flavor unpolarized par- 
ton distributions. It was found that the predicted average transverse momentum square of 
quarks and antiquarks depends strongly on their longitudinal momentum fraction x. We 
also estimate the average transverse momentum of quarks and antiquarks separately to find, 
somewhat unexpectedly, that (kl)^ ~ 0.224 GeV^ and {kl)^ ~ 0.445 GeV^ that is, the 
antiquarks have much higher average transverse momentum than the quarks. On the other 
hand, the average momentum transfer of quarks and antiquarks altogether turns out to be 
{kj_)^~^^ ~ 0.266 GeV^, which is order of magnitude consistent with the recent phenomeno- 
logical analysis, although we must be careful about the fact that the average transverse 
momentum square is a Q^-dependent quantity. 
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